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The cross section for e + e~ — > rjJ/ip between yfs = 3.8 GeV and 5.3 GeV is measured via initial 
state radiation using 980 fb -1 of data on and around the T(nS)(n = 1, 2, 3, 4, 5) resonances collected 
with the Belle detector at KEKB. Two resonant structures at the '0(4040) and t/>(4160) are observed 
in the r\J '/ip invariant mass distribution. Fitting the mass spectrum with the coherent sum of 
two Breit-Wigner functions, one obtains 23(i/>(4040) -> rjJ/ip) ■ T^ 040 ' = (4.8 ± 0.9 ± 1.4) eV 
and 8(V>(4160) -»■ rjJ/tp) ■ I^ 4160 ' = (4.0 ± 0.8 ± 1.4) eV for one solution and 23(V>(4040) -»■ 

rjJ/ip) ■ r^ 4040 > = (11.2 ± 1.3 ±1.9) eV and B (if) (4160) -> rjJ/ip) ■ r^ 41 _ 60) = (13.8 ± 1.3 ± 2.0) eV 
for the other solution, where the first errors are statistical and the second are systematic. This is 
the first measurement of this hadronic transition mode of these two states, and the partial widths 
to r)J/if> are found to be about 1 MeV. There is no evidence for the V(4260), y(4360), if) (4415), or 
y(4660) in the T]J/ip final state, and upper limits of their production rates in e + e" annihilation are 
determined. 

PACS numbers: 14.40.Pq, 13.25.Gv, 13.66. Be 



Many charmonium and charmoniumlike states have 
been discovered at B-factories in the past decade. Some 
of these states are good candidates for conventional char- 
monium states, while others exhibit unusual properties 
consistent with expectations for exotic states such as a 
multi-quark state, molecule, hybrid, or the glueball [lj. 
In the vector sector, four exotic charmoniumlike struc- 
tures, F(4008) and V(4260) in e+e" -> tt+tt" J/ip BM 
and F(4360) and V(4660) in e+e" -> tt+tt-^S) fflE 
have been reported via initial state radiation (ISR), in 
addition to the three known excited ip states above 
4.0 GcV/c 2 : V(4040), V(4160), and ^(4415). It is un- 
likely that all seven of these states are charmonia, as 
the potential models predict only five vector states in 
this mass region . The current understanding of these 
states is based on limited statistics, and the fact that 
some may be produced via mechanisms that are diffi- 
cult to estimate theoretically, such as final state rescat- 
tering [lj, makes the determination of which might be 
exotic even more challenging. In order to further the un- 
derstanding of the nature of these states, it is important 
to investigate them using much larger data samples. 

An important study is the investigation of hadronic 
transitions (either by an r/ or a pion pair) between 
these states and a lower charmonium state like the J /ip. 
The CLEO collaboration measured er(e + e~ — » rjJ/tp) = 



I5±l ± 8 pb at y/s = 4120 - 4200 MeV Q, and the BE- 
SIII collaboration reported cr(e + e~ — > rjJ/ip) — (32.1 ± 
2.8) pb at y/s = 4009 MeV [8J, which is in agreement 
with the theoretical calculation including contributions 
from the known if) states and the virtual charmed meson 
loops Q. However, the limited statistics of the CLEO 
analysis prevented the measurement of the line shape 
of r]J/ip. Thus, it is worthwhile to study the process 
e + e~ — > r/J/tp via ISR with the full Belle data sample to 
search for rj transitions from these seven states to J ' /if). It 
is worth noting that the ip states arc identified in decays 
to charmed meson pairs but not in dipion transitions to 
lower if; states, while the opposite is true of the Y states. 
There may also be surprises from transitions of unex- 
pected states. 

In this Letter, we report an investigation of the 
e + e~ — > r/J/ip process using ISR events observed with 
the Belle detector [l(| at the KEKB asymmetric-energy 
e + e~ collider [ll|. Here, J/ip is reconstructed in the 
£ + £~ (I = e, /i) final state and -q in the 77 and tt + tt^tt° 
final states. Due to the high background level from 
Bhabha scattering, the J/ip — > e + e~ mode is not used 
in conjunction with the decay mode rj — > 77. The inte- 
grated luminosity used in this analysis is 980 fb _1 . About 
70% of the data were collected at the T(45) resonance, 
and the rest were taken at other T(nS) (n = 1, 2, 3, or 
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5) states or center-of-mass (CM) energies a few tens of 
McV lower than the T(45) or the T(raS) peaks. 

We use the PHOKHARA event generator [12j to simulate 
the process e + e _ — > ^isRjjJ/ip. In the generator, one 
or two ISR photons may be emitted before forming the 
resonance A, which then decays to r]J/ip, with J/ip — > 
e + e~ or and 77 — > 7r + 7r~7r° or 77. 

For a candidate event, we require two (four) good 
charged tracks with zero net charge for 77 — > 77 (rj — > 
7r + 7T~7r°). A good charged track has impact parameters 
with respect to the interaction point of dr < 0.5 cm in 
the r-(j> plane and \dz\ < 5 cm in the r-z plane. The 
transverse momentum of the leptons is required to be 
greater than 0.1 GeV/c. For each charged track, infor- 
mation from different detector subsystems is combined 
to form a likelihood for each particle species (i), Li [Hj]. 
Tracks with TZk = cT+c < ^ are identified as pi- 
ons with an efficiency of about 95%, while 6% of kaons 
are misidentified as pions. Similar likelihood ratios are 
formed for electron and muon identification 3, IH ■ For 
electrons from J/ip — > e + e~ , both tracks are required to 
have 7Z e > 0.1. The bremsstrahlung photons detected in 
the electromagnetic calorimeter (ECL) within 0.05 radi- 
ans of the original e + or e~ direction are included in the 
calculation of the e + e~(7) invariant mass. For muons 
from J/ip —> one of the tracks is required to have 

TZfi > 0.9 and the other track should have associated hits 
in the if z,-and-muon detector (KLM) that agree with the 
extrapolated trajectory of a charged track provided by 
the drift chamber. The lepton ID efficiency is about 90% 
for J/ip -> e+e - and 87% for J/ip -> 

The 77 is reconstructed from 7r + 7r~7r and 77 final 
states. To reconstruct 77 — > 7r + 7r~7r°, the ir° is recon- 
structed from two photons. A photon candidate is an 
ECL cluster with energy ^(7) > 25 MeV that does not 
match any charged tracks. The mass resolution of ir° is 
about 5.2 McV/c 2 from MC simulation. Considering the 
low-mass tail, the invariant mass of the photon pair is 
required to be between 110 MeV/c 2 and 150 MeV/c 2 for 



a 7r candidate. 



combinations are formed and 



are subject to a mass-constrained kinematic fit. When 
there is more than one tt° candidate, the combination 
with the smallest % 2 from the mass-constrained fit is se- 
lected as the 77 candidate. Events with 7-conversions are 
removed by requiring lZ e < 0.75 for the it + it~ tracks 
from 77 decays. In the reconstruction of 77 — > 77 candi- 
dates, two photon candidates are required with energies 
in the laboratory frame satisfying E(-fi) > 0.15 GeV and 
E(jh) > 0.4 GeV, where the subscript I (h) signifies the 
lower (higher) energy photon. 

The scatter plots of dilepton invariant mass M^+i- 
versus 77-candidate invariant mass M„+ w - n o or 7/7^ 
invariant mass Af 77 are shown in Fig. [1] for events 
that survive these selection criteria. Here the invari- 
ant masses are calculated with the momenta before the 
mass constraints. A dilepton pair is considered as a 



J/ip candidate if M e+e - is within ±45 MeV/c 2 (the 
mass resolution being 15 MeV/c 2 ) of the J/ip nominal 
mass. The J/ip mass sidebands are defined as Mg+£- G 
[3.172,3.262] GcV/c 2 or M i+i - G [2.932,3.022] GcV/c 2 . 
A fit of the M^+ w - n o or M 77 distribution with a Gaus- 
sian plus a second-order polynomial yields a mass res- 
olution of 4.3 MeV/c 2 for the 77 — > ir + Tr~ir° mode and 
11.1 MeV/c 2 for the 77 — > 77 mode. We define the 77 
signal region as M^+^-^a G [0.5343,0.5613] GeV/c 2 and 
M 77 G [0.5, 0.6] GeV/c 2 and the 77 mass sideband regions 
as M„+ w - v o G [0.5748,0.6018] GcV/c 2 or M w+v - n o € 
[0.4938,0.5208] GeV/c 2 , and M 77 G [0.35,0.45] GcV/c 2 
or M 77 G [0.65,0.75] GcV/c 2 . The central (surround- 
ing) rectangles of Fig. [1] show the rjJ/ip signal (side- 
band) regions. With 51 (52) representing the sum of 
the events in the four sideband boxes nearest (diagonal) 
to the signal box, the normalization of the sidebands is 
5 = 0.5 x 51 - 0.25 x 52. 
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FIG. 1: Invariant mass distributions of (a) £ + £~ vs. 7r + 7r~7r° 
and (b) £ + £~ vs. 77 for selected ir + ir~ir °£ + £- or 77^+^- 
candidates with invariant mass between 3.8 GeV/c 2 and 
5.3 GeV/c 2 . The box in the center of each plot shows the 
r/J/tp signal region while the surrounding boxes show the side- 
band regions. 



The detection of the ISR photon is not required; 
instead, we require —1 (GcV/c 2 ) 2 < M 2 ec < 
2.0 (GeV/c 2 ) 2 , where M 2 ec is the square of the mass re- 
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coiling against the rjj/i/j system. In calculating M r 2 0C , 
the momenta of the J/ip and 77 after the kinematic fit are 
used to improve the resolution of Af 2 ec . The fit constrains 
signal candidates to the 77 and J/ip masses, while events 
having r\ or J/ip candidate masses lying in sideband re- 
gions are fitted with masses constrained to the center of 
the sideband region. 

Figure [2] shows the 77 J/ip invariant mass (M vJ ^ (lij ) 
for selected candidate events, together with background 
estimated from the scaled 77 or J/ip mass sidebands. Two 
distinct peaks are evident in Fig. f2J one at 4.0 GeV/c 2 
and the other at 4.2 GcV/c 2 , in addition to the dominant 
ip(2S) signal. The cross section of e + e~ — > 7ISR'0(2»S , ) 
in the full Belle data sample is measured to be 13.9 ± 
1.4 (stat.) pb using the 186 ± 17 77 —> 7T + 7r~7r° events 
and 14.0 ± 0.8 (stat.) pb using the 470 ± 25 77 -> 77 
events, in good agreement with the production cross sec- 
tion of 14.7 pb calculated by using the world average 
values of the mass, width, and partial width to e + e~ of 
ip(2S) [l7| . and the e + e~ CM energies correspond to the 
Belle data samples. Agreement between the J/ip — > e + e~ 
and J/ip — > modes is also observed. The vis- 

ible energy (_E v is) an d polar angle distributions of the 
r)J/tp system in the e + e~ CM frame for the events with 
M nJ /^ G [3.8,5.3] GcV/c 2 agree well with the MC sim- 
ulation, confirming that the signal events are produced 
via ISR. Here, E v [ s is the total energy of all final state 
photons and charged particles. Charged particle energies 
are calculated from track momenta assuming the tracks 
to be pions. 

An unbinned maximum likelihood fit is performed to 
the mass spectra M v j/^ G [3.8,4.8] GeV/c 2 from the 
signal candidate events and 77 and J/ip sideband events 
simultaneously, as shown in Fig. [3J The fit to the sig- 
nal events includes two coherent P-wavc Breit-Wigncr 
functions, BW 1 for 0(4040) and BW 2 for 0(4160), as- 
suming that only two resonances contribute to the nJ/tp 
final states, and an incoherent second-order polynomial 
background; the fit to the sideband events includes the 
same background function only. The width of each res- 
onance is assumed to be constant, and an overall two- 
body phase-space factor is applied in the partial width 
to rjJ/tp. The signal amplitude is M = BWi + ■ BW 2 , 
where <p is the the relative phase between the two reso- 
nances. In the fit, the BW functions are convolved with 
the effective luminosity [18[ and M^j/^-dependent effi- 
ciency, which increases from 4% at M^j/^ = 4.0 GeV/c 2 
to 7% at M v j/,p = 4.5 GeV/c 2 . The effect of mass res- 
olution, which is determined from MC simulation to be 
5 — 11 MeV/c 2 over the resonant mass region, is small 
compared with the widths of the observed structures, 
and therefore is neglected. A fit performed with floating 
masses and widths for the two structures yields a mass 
of (4012 ± 5) MeV/c 2 and width of (54 ± 13) MeV for 
the first, and a mass of (4157± 10) MeV/c 2 and width of 
(84 ± 20) MeV for the second. Their masses and widths 
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FIG. 2: The invariant mass distribution of the rjJ/tp candi- 
dates. The top row shows the rj — > tt + tt~tt° mode and the 
bottom row shows the rj — > 77 mode. The open histograms 
are from the 77 and J/ip signal region, while the shaded ones 
are from their sideband regions after the proper normaliza- 
tion. The insets show the distributions around the ip(2S) 
mass region. 



arc in agreement with those of the -0(4040) and 0(4160), 
and thus they are referred to hereafter as the -0(4040) 
and -0(4160). In the fit below, the masses and widths 
of these two resonances are fixed to their world average 
values [TtJ as the statistics are low here. 

Figure[3]and Tablcfl] 19] show the fit results. There are 
two solutions with equally good fit quality. To determine 
the goodness of the fit, we bin the data (events in both 
signal and sideband regions) so that the expected number 
of events in a bin is at least seven and then calculate a 
X 2 /ndf of 71.4/46, corresponding to a confidence level 
(C.L.) of 0.9%, where ndf is the number of degrees of 
freedom. The significance of each resonance is estimated 
by comparing the likelihood of fits with and without that 
resonance included. We obtain a statistical significance 
of 6.5cr for -0(4040) and 7.6ct for ?0(416O). Varying the 
masses and widths of resonances by lcr, the fit range by 
200 MeV/c 2 , and the order of the background polynomial 
by one, we obtain a minimum statistical significance of 



5 



6.0cr for ^(4040) and 6.5er for ^(4160). 

Taking r^ 040) = (0.86±0.07) keV from PDG 0, one 
obtains 23(0(4040) -> r]J/tp) = (0.56 ± 0.10 ± 0.17)% or 
23(0(4040) -» r)J/ip) = (1.30 ±0.15 ±0.24)%; while using 
the PDG average value r^f_ 60) = (0.83 ±0.07) keV 0, 
one gets 23(0(4160) r?J/^) = (0.48 ± 0.10 ± 0.17)% 
or (1.66 ± 0.16 ± 0.28)%. In each case, the first error is 
statistical and the second is systematic. These indicate 
the transition rates of these states to r]J/ip are large, 
being of order 1 McV. 

Possible contributions from other excited charmo- 
nium(likc) states are examined. There is a cluster of 
events near the M n ji^ = 4.36 GcV/c 2 . Assuming it is 
the F(4360), the significance is l.ler in a fit with the 
masses and widths of the 0(4040) and ■0(4160) fixed to 
their world average values [l7j, or 2.9ct if the masses 
and widths of 0(4040) and 0(4160) are free. Besides 
the F(4360), the F(4260), 0(4415) and y(4660) are in 
[3.8,5.3] GeV/c 2 mass region. Fits that include each 
one of them and the masses and widths of 0(4040) and 
0(4160) fixed to their world average values [17J are per- 
formed to determine the upper limits of 23 ■ T e + e - . The 
systematic errors that will be described later in the text 
together with those from the uncertainties of the 0(4040) 
and ■0(4160) resonant parameters are considered in the 
upper limit determination. In order to be conserva- 
tive, the efficiencies have been lowered by a factor of 
1 — a sys in the calculation. We obtain the upper lim- 
its on B(X -> ? ? J/0) • Tf +e _ for X = F(4260), y(4360), 
0(4415) and F(4660) are 14.2 eV, 6.8 eV, 3.6 eV and 
0.94 eV at 90% C.L., respectively. 




M(r|J/\j/) (GeV/cf) 



FIG. 3: The r/J/if) invariant mass distribution and the fit re- 
sults. The points with error bars show the data while the 
shaded histogram is the normalized r\ and J/0 background 
from the sidebands. The curves show the best fit on signal 
candidate events and sideband events simultaneously and the 
contribution from each Breit-Wigner component. The inter- 
ference between the two resonances is not shown. The dashed 
curves at each peak show the two solutions (see text). 



TABLE I: Results of the fits to the ??J/0 invariant mass spec- 
trum. The first errors are statistical and the second are sys- 
tematic. M, F, and B ■ r* +e _ are the mass (in MeV/c 2 ), 
total width (in MeV), product of the branching fraction of 
V' — > riJ /ijj and the — > e + e~ partial width (in eV), respec- 
tively. is the relative phase between the two resonances (in 
degrees). 



Parameters Solution I 



Solution II 



M, 



^(4040) 
^(4040) 



B-V 



i/>(4040) 



M, 



^(4160) 



- V(4160) 



B-V 



i/>(4040) 



4039 (fixed) 
80 (fixed) 
4.8 ±0.9 ±1.4 11.2 ±1.3 ±1.9 
4153 (fixed) 
103 (fixed) 
4.0 ±0.8 ±1.4 13.8 ±1.3 ±2.0 
336 ±12 ± 14 251 ±4±7 



To estimate the errors in B ■ T e + e ~ , the uncertainties 
from the choice of parametrization of the resonances (es- 
pecially introducing the mass dependence for the widths), 
the masses and widths of resonances 17J , the fit range, 
the background shape and the possible contributions 
from 0(25) or 0(4415) are considered. The total errors 
are 35.0% and 14.8% for solutions I and II, respectively. 
The particle ID uncertainty is 5.5%; the uncertainty in 
the tracking efficiency is 0.35% per track and is additive; 
the uncertainty in the photon reconstruction is 2% per 
photon. The uncertainties in the J/0 mass, r\ mass, and 
A/ t 2 cc requirements are measured with the control sample 
e + e~ — > 0(25) — > ?yJ/0. The efficiencies of the require- 
ments on the data are obtained from the fits of the corre- 
sponding distributions. The MC efficiency is found to be 
higher than in data by (2.3±2.6)% for the 7r + 7r~7r° mode 
and (0.1 ± 1.6)% for the 77 mode. A correction factor 
1.023 is applied to the 7r + 7r~7r° final state, and 2.6% is 
conservatively taken as the associated systematic error of 
the sum for 7r + 7r~7r° and 77 modes. 

Belle measures luminosity with 1.4% precision while 
the uncertainty of the generator phokhara is less than 



1% [I2J. The trigger efficiency for the events surviving 



the selection criteria is around 91% with an uncertainty 
smaller than 2%. The uncertainties in the intermediate 



decay branching fractions taken from Ref. [17[ contribute 
a systematic error of less than 1.6%. The statistical error 
in the MC determination of the efficiency is 0.2%. 

Assuming all the sources are independent and adding 
them in quadrature, we obtain total systematic errors in 
B ■ r e+e - of 36% for Solution I and 17% for Solution II 
for both V(4040) and ^(4160). 

The cross section for e + e~ — > r/J/ip for each r\J /ip mass 
bin is calculated according to 
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where j is the j-th mode of r/J '/ip decays (j 



"7r°e+e- 



~ir 7r°/i + /i , and 77/i + /^ ); n[ 



obs n bk s 



£y, and Bj are number of events observed in data, 
number of background events estimated from sidebands, 
detection efficiency of the j-th mode, effective luminosity 
in the i-th r]J/tjj mass bin, and the branching fraction of 
r/J/ip decays into the j-th mode [T3|, respectively. The 
resulting cross sections in the full solid angle arc shown 
in Fig. IH where the error bars include the statistical un- 
certainties in the signal and the background subtraction. 
The systematic error for the cross section measurement, 
which includes all the sources that have been described 
other than those arising from the details of the fit to 
the mass spectrum, is 8.0% and common to all the data 
points. The cross sections of e + e~ -4 i]J/tp are around 
70 pb and 50 pb at the ■0(4040) and ^(4160) peaks, re- 
spectively, to be compared with around 20 pb and 10 pb 
measured in e + e~ — > tt + it~ J/xJj 




FIG. 4: The measured e + e 



r\J jip cross section for ^/s ■■ 



3.8 GeV to 5.3 GeV. The errors are the summed statistical 
errors of the numbers of signal and background events. A 
systematic error of 8.0% common to all the data points are 
not shown. 

In summary, the e + e~ — > rjj/^fj cross section is mea- 
sured from 3.8 GeV up to 5.3 GeV for the first time. Two 
distinct resonant structures, the -0(4040) and 0(4160), 
are observed. This is the first time that the 0(4040) 
and '0(4160) have been observed to decay to final states 
not involving charm meson pairs. The products of the 
branching fraction to rjj/ip and the e + e~ partial width 

V>(4040) 



are determined to be S(^(4040) -> i]J/ip) ■ T 

(4.8 ± 0.9 ± 1.4) eV and 6(0(4160) -> r]J/ip) ■ F 
(4.0 ± 0.8 ± 1.4) eV for one solution; or £(■0(4040) 



</>(4160) 



riJH>) ■ 

r,J/i>) 



^(4040) 
^(4160) 



(11.2±1.3±1.9)eV and B(V>(4160) 
= (13.8 ± 1.3 ± 2.0) eV for the other 



solution. These transition rates correspond to about 
1 MeV partial widths to rjj/ip for these two states. We 
find no evidence for the V(4260), V(4360), V(4415) or 



V(4660) in the i]J/ip final states, and upper limits of 
their production rates in e + e~ annihilation are deter- 
mined. The present measurement reveals clear peaks due 
to the 1^(4040) and 0)(416O) decays observed in experi- 
mental data that are absent in the prediction in Ref. 0] , 
although the theoretical calculation with carefully cho- 
sen parameters agrees with the measured cross sections 



of e H 



r/J/tp. 
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